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Abstract. Plateau de Bure Interferometer high-spatial resolution CO observations combined with near-infrared 
H2 data disentangle at least seven (maybe even nine) molecular outflows in the massive star-forming region 
IRAS 19410+2336. Position-velocity diagrams of the outflows reveal Hubble-like relationships similar to outflows 
driven by low-mass objects. Estimated accretion rates are of the order 10~ 4 Mq yr" 1 , sufficiently high to overcome 
the radiation pressure and form massive stars via disk-mediated accretion processes. The single-dish large-scale 
mm continuum cores fragment into several compact condensations at the higher spatial resolution of the PdBI 
which is expected due to the clustering in massive star formation. While single-dish data give a simplified picture 
of the source, sufficiently high spatial resolution resolves the structures into outflows resembling those of low- 
mass star-forming cores. We interpret this as further support for the hypothesis that massive stars do form via 
disk-accretion processes similar to low-mass stars. 
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1. Introduction 

The physical processes forming massive stars are still sub- 
ject to great debate. Basic low-mass star-forming theo- 
ries predict small accr etion rates independent of mass 
(lO^ 6 - lO- 5 M yr-\ Ishul Il977|b which are not suffi- 
cient to overcome the radiation pressure of a star > 
8M Q . More massive stars should not fo rm according to 
that picture l|Wolfire fc Cassinellil Il987ft although they 
are observed to exist. Several theories have been put 
forward to solve this discrepancy: if the standard sce- 
nario is adapted, higher accretion rates must be present 
and the accretion is likely to be mediated through 



disks 


(e.g., 


Wolfire k Cassinellil Il987t Ijiiina fc Adams 


1996; 


Norberg & Maeden 


2000; 


iMaeder & Behrendl 


l2002t 
lother 


Yorke & Sonnhalterll2002tlMcKee & TarJl2002t. Ai 



scenario suggests that the very center of the evolving clus- 
ter is dense enough that intermediate-mass protostars col- 
lide, merge, and form the most mas s ive objects via these 
collisions (e.g.. iBonnell et alJll998t IStahler et all l2000t 
Ballvl l200a IZinnecker fc BatdboO: 
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As high-mass star-forming regions are on the average 
far away (a few kpc) and massive stars form in a clus- 
tered mode, it is nearly impossible to resolve with cur- 
rent telescopes the forming cluster in the mm and sub- mm 
regime, and thus to study the physical processes forming 
the stars in a direct way. Infrared observations give am- 
biguous results at best, since the centers of activity are 
very deeply embedded and block out all infrared light, ex- 
cept in rare favorable geometries. Therefore, indirect ap- 
proaches have to be taken. Massive molecular outflows are 
perfectly suited for such investigations, because outflows 
are observed on large spatial scales (in the parsec regime), 
they are easier to resolve spatially, and their structure and 
energetics give information on the physical processes close 
to their driving center. For example, the accretion scenario 
predicts outflows which are morphologically similar to low- 
mass outflows and which show a high degree of collimation 
due to the star-disk interaction ijRicher et all 12000). On 
the contrary, colliding protostars are expected to be ex- 
tremely eruptive phenomena during which accretion disks 
should not be able to survive, and therefore the outflows 
are likely less collimated and rather resemble explosions 
as observed in Orion KL jSchultz et allll999|b 

Over recent years several single-dish studies 
of massive outflows have b een undertaken (e.g., 
Shepherd fc Churchwelll Il996l iRidge fc Moorel l200ll 



Zha ng et all l200lt iBeuther et all l2002b(l . While these 



studies agree on the facts that outflows are ubiquitous 
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phenomena in massive star formation and that they are 
very massive and energetic, there is considerable disagree- 
ment on the typical coll im ation of the ob served outflows. 
ISheoherd fc Churchwelll il99(j\ and iR.idge fc MooTel 
l)200l[) claim that the average collimation of massive 
outflows is lower than observed for their low-mass 
counterparts supporting the proposal that massive star 
formation has to proceed through diff erent physica l 
processes such as coalescence. iBeuther et, all ( 20f)2h|) 
also observed lower collimation of massive outflows 
compared to low-mass sources. However, taking properly 
into account the poor spatial resolution, they argue that 
their data are well consistent with highly collimated 
outflows in high-mass star-forming regions (but proof of 
this requires interferometer observations). Thus, no other 
physical processes have to be invoked, and high-mass 
star formation can proceed as in the classical low-mass 
scenario, just with significantly enhanced accretion rates. 

Despite the great advances during recent years, single- 
dish observations are not sufficient to understand the 
complex outflows in massive star formation. As always 
clusters of stars form simultaneously, it is likely that 
at high spatial resolution many of the single-dish out- 
flows resolve into multiple outflows emanating from the 
same large scale core but from different protostars within 

them. Using P l ateau d e Bure interferometric observations 

| n I i — i 
Beu ther et alJ il2002di showed that the chaotic single-dish 

outflow in IRAS 05358+3543 resolves into at least three 
molecular outflows at a spatial resolution of ~ 3.5". One of 
these outflows shows a collimation degree of 10, the high- 
est so far observed in massive star-forming regions, and 
as high as the highest values known for low-mass star for- 
mation. The observations indicate that high-mass regions 
are more complex due to the clustered mode of formation, 
but that the physical processes taking place are similar to 
those of their low-mass counterparts. 

There do exist other interferometric observa- 
tions of high-mass o utflows, e.g., IRAS 2 0126+2104 
(ICe^aroifi^t^alJ, |l997j I S hepher d et_al 1 12000(1. G1 92.16 
llSheoherd et*all Il998t) . G35.2 llOibb et all l2003h . but 



they are either covering only the very inner region not 
tracing the large scale outflow, or the spatial resolution 
does not exceed 5". 

In order to base the results obtained so far on a 
larger statistical base and observe more massive out- 
flows at adequate high spatial resolution, we perform a 
series of high-resolution observations using the Plateau 
de Bure Interferometer (PdBI) and BIMA. Here, we 
present the results of PdBI observations of the outflow 
system in IRAS 19410+2336 which is part of a sam- 
ple of 69 high-mass protostellar objects (HMPOs) stud- 
ied oyer recent years in grea t detail I Sridharan et all 
l2002t IBeuther et all l2002albllcl l. At the very center of 
the southern core, we detect H2O and CH3OH maser 
emission as well as a weak 1 mJy unresolved cm con- 
tinuum source at a spatial resolution of 0.7" (see Fig. 

Assuming the latter to be due to optically thin free- 
free emission at 3.6 cm, it indicates a recently ignited 



but not very evolv e d mas s ive object at the clust er cen- 
ter (IBeuther et all l2002ct ISridharan et all l2002|) . It is 
also the most deeply embedded massive star-forming re- 
gion whi ch was ever observed w ith the CHANDRA X-ray 
satellite ( Beuther et all l2002ejl . So far, we only knew the 
kinemat ic distances which suffer from a distance ambigu- 
ity l|Sridharan et all I2OO2]) . but based on the new MSX 
and 2MASS data of the region we believe the near kine- 
matic distance to be more likely (Bontemps, priv. comm.). 
Further on, we assume IRAS 19410+2336 to be at a dis- 
tance of ~ 2.1 kpc. The approximate bolometric luminos- 
ity of the region is 10 4 L Q . Single-dish 1.2 mm continuum 
observation reveal two massive adjacent dust cores, each 
associated with a bipolar molecular outflo w as observed 
in CO (2 -1) with 11" resolution (Fig. □ IBeuther et all 
2002a b). The single-dish data suggest a simple morphol- 
ogy in both cores which makes the region a good candidate 
for further studies of collimated outflows. The PdBI obser- 
vations will show that the single-dish observations provide 
only a simplified picture of the real outflow structure. 



2. Observations 

2.1. Plateau de Bure Interferometer (PdBI) 

We observed IRAS 19410+2336 in Summer 2001 with 
the Plateau de Bure Interferometer at 2.6 mm in the 
D (wit h 4 antennas) and C ( with 5 antennas) configu- 
ration l)Guilloteau et all Il992h . The simultaneously ob- 
served 1 mm data were only used for phase corrections be- 
cause of the poor Summer weather conditions. The 3 mm 
receivers were tuned to 115.27 GHz (USB) (centered at 
the 12 CO 1 — > line) with a sideband rejection of about 
5dB. At this frequency, the typical SSB system tempera- 
ture is 300 to 400 K, and the phase noise was below 30°. 
The frequency resolution at 115GHz was O.lkms -1 , but 
we smoothed the data to lkms -1 which is sufficient to 
map the broad line wings. Atmospheric phase correction 
based on the 1.3 mm total power was applied. For con- 
tinuum measurements, we placed two 320 MHz correla- 
tor units in the band to cover the largest possible band- 
widths. Temporal fluctuations of the amplitude and phase 
were calibrated with frequent observations of the quasars 
1923+210 and 2023+336. The amplitude scale was derived 
from measurements of MWC349 and CRL618, and we es- 
timate the final flux density accuracy to be ~ 15%. To 
cover both cores a mosaic of 10 fields was observed. We 
obtain a 3er continuum rms o f ~ 5 mJy. The contin uum 
data are already presented bv IBeuther et alJ 1 2002eh . 



2.2. Short spacings with the IRAM 30m telescope 

To account for the missing short spacings and to recover 
the line-flux, we also observed the source in CO (1-0) with 
the IRAM 30 m telescope in Summer 2002. The obser- 
vations were done remotely from Bonn in the on-the-fly 
mode. Due to previous software changes the online veloc- 
ity correction was missing a term corresponding to some 
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part of the earth rotation. This error was corrected offline 
by the AoD Gabriel Paubert. The velocity is estimated to 
be correct within 0.1 km s -1 . 

The algorithm to derive visibilities from the single-dish 
da ta corresponding t o each pointing center is described 
bv lGueth et alllll996h . The single-dish and interferometer 
visibilities are subsequently processed together. Relative 
weighting has been chosen to minimize the negative side- 
lobes in the resulting dirty beam while keeping the high- 
est angular resolution possible. Images were produced us- 
ing natural weighting, then a CLEAN-based deconvolu- 
tion of the mosaic was performed. The final beam size of 
the 2.6mm data is 3.9" x 3.6" (P.A. 33°). 

2.3. CalarAlto 

The near-infrared camera Omega Prime on the 3.5 m tele- 
scope on Calar Alto/Spain was used to obtain 2.12/im 
narrow band H2 and K' wide field images of the region 
around IRAS 19410+2336 on June 11 2001. At a pixel 
scale of 0.4" pixel -1 , the 1024x1024 pixel array provides 
a field-of-view of ~ 6"8 x 6" 8. The basic data redu ction 
has already been described in iBeuther et all l)2002e|) . The 
narrow-band line mosaic and the K' mosaic were then 
carefully registered. The Point Spread Function (PSF) of 
the H 2 image was degraded to match the PSF of the K' 
image with the poorer seeing; finally, the continuum im- 
age was scaled such that stars had the same flux as in 
the narrow band image. This image was then subtracted 
from the (PSF-degraded) narrow band image in order to 
provide the continuum subtracted narrow band image (see 
Figure |3 for the K' and H 2 image). 

3. Results 

3.1. Millimeter continuum sources 

The large-scale southern single-dish core splits up into at 
least four sub-cores in the 2.6 mm continuum (Fig. QJ, 
whereas the northern core is resolved into 2 sub-sources in 
the higher-resolution data. Assuming that the 2.6 mm con- 
tinuum is produced by optically thin thermal dust emis- 
sion, we can calculate the masses a nd the column dens ities 
of the different sub-cores following iHildebrandl (^983): 



1.3 x 10 -3 a p R F v 
JJT) 0.1^m3gcm -3 100 Jy 




7.8 x 10 1 " a p R F v 
J v (T)n 0.lAim3gcm- 3 100 Jy 




where J V {T) = [cxp(hv/kT) — l] -1 and f2,a,p, R, and (3 
are the beam solid angle, grain size, grain mass density, 
gas-to-dust ratio, and grain emissivity index for which we 
use the values 3.9" x 3.6" in radians, 0.1 /im, 3gcm~ 3 , 100, 



and 2, respectively l)Hunter et allEoOOl) . Based on IRAS 
far-infrared data, the dust tempe rature T is assumed t o 
be around 30 K. As discussed bv IBeuther et alJ l(2002ah . 
higher dust opacity indices (corresponding to a lower (3) 
would result in lower masses and column densities whereas 
lower temperatures would increase the derived parame- 
ters. We estimate the results to be correct within a factor 
5. 




19 h 43 m 10 s 19 h 43 m 12 s 10 : 

R.A. [J2000.0] 



Fig. 1. The left panel p resents the single-dish observa- 
tions already published in lBeuther et alJ l|2002tJ) . In grey- 
scale we show the 1.2 mm dust emission, and the solid 
and dashed contours outline the blue and red CO(2-l) 
outflow emission. Both datasets were observed with the 
IRAM 30 m telescope at a spatial resolution of 11" (bot- 
tom left). The arrows mark the two outflows identified 
previously from these observations. In the right panel, 
we present the PdBI 2.6 mm continuum data and label 
the sub-sources. The beam is presented at the bottom 
right and the white cross marks the position of the cm 
continuum source. 

The masses range between 40 and 180 M© and the col- 
umn densities are in the 10 24 cm -3 regime, corresponding 
to a visual extinction of A v ~ 1000 (see Table l3~T)l . The 
mm sources mml and mm2 are of similar mass but the 
most prominent sub-core is mml, which coincides with 
the H 2 and CH3 OH masers and also w ith the weak cm 
continuum source ijBeuther et al.ll2002ej) . The continuum 
3a rms is ~5mJy equal to a mass sensitivity >20M© at 
a temperature of 30 K. Assuming an S v ex v A relationship 
(i.e., (3 — 2) for the continuum dust emission in the mm 
regime, we compare the single-dish flux at 1.2 mm with the 
PdBI data at 2.6 mm and find that approximately 60% 
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of the total continuum flux is filtered out by the inter- 
ferometric observations. This most likely corresponds to 
emission from the large-scale envelope of the core, whereas 
the PdBI detected mm emission traces the very dense and 
compact condensations. As the weak 1 mJy cm continuum 
emission is most likely due to optically thin free- free emis- 
sion, its contribution in the mm-regime is neglectable. 

3.2. Outflow morphology 

Figure [3 shows a channel map of the merged PdBI+30m 
CO emission in IRAS 19410+2336, and we find a multi- 
tude of different features. Based on the channel maps the 
velocity-regimes of outflowing gas are defined as: blue wing 
[5;18] km/s and red wing [26;47] km/s. The velocity spread 
does not vary significantly for the different sub-outflows 
described below. 

Comparing the single-dish CO emission presented in 
Figure with the high-resolution outflow map of the 
PdBI (Fig. |3| (a) and (b)) we see that the overall struc- 
tures in both datasets are similar, we just observe - as 
expected - much more structure in the interferometric 
data. Identifying the same structures makes us confident 
that the data reduction and merging process has been 
done properly, but the sub-structures demand a different 
morphological interpretation of the molecular outflows. 
Additionally, the near-infrared shocked H2 data show 
many emission knots associated with the star-forming re- 
gion (Fig. Most of the H2 features can be correlated 
with the molecular gas, but for some H 2 knots this is less 
obvious. 

We identify four molecular outflows in the southern 
core and three outflows associated with the northern core. 
The identifications of the outflows are not unambiguous, 
and we are still resolution limited. Therefore, we note that 
the identifications below are just one interpretation, oth- 
ers might be possible and only very-high-spatial-resolution 
observations anticipated for next winter will clarify the 
picture. 

3.2.1. The southern region 

The morphological interpretation of the four southern out- 
flows is sketched in Figures [21 and [31 

South (A): The structures labeled (A) and (B) could 
also be part of only one outflow, but combining the CO 
and H2 observations this does not seem likely and we pre- 
fer the interpretation of two different overlapping outflows. 
The outflow South (A) corresponds to the main large scale 
outflow already identified in the single-dish observations 
(see Figures ^ 121 & ED ■ It runs in east- west direction and 
is centered approximately on the main southern sub-core 
mini which also houses the cm-source and the maser emis- 
sion. Morphologically, one of the X-ray sources could also 
be the center of the outflow, but the main mm source 
seems to be more likely. 



South (B): The H 2 features (1, 4, 5, 6) and the cor- 
responding CO emission outline the ellipsoidal structure 
of outflow South (B) sketched in Figures [21 and [3J As the 
molecular emission at the eastern end of the outflow has a 
strong gradient and stops rather instantaneously directly 
before the H2 feature (1), it is likely that both - CO and 
H2 - outline a bow-shock-like working surface between an 
underlying jet and the surrounding molecular gas. The 
mm source mm3 likely houses the driving source of out- 
flow South (B). As discussed below, the H 2 features (4), 
(5) and (6) might outline the interface between the outflow 
cavities of outflow (B) and the northern outflow (F). 

South (C): The third southern outflow (C) is virtu- 
ally undetectable in the single-dish data because it is very 
small in extent but clearly resolved in east-west direction 
in the high-resolution PdBI data. The elongated H 2 emis- 
sion feature (2) is strongly associated with the blue molec- 
ular gas and outlines the same morphological structure, 
indicating a highly collimated outflow. The H 2 knot (3) 
is at the eastern end of the red CO gas and can be part 
of this outflow as well. South (C) is near the center of 
the region and might also be driven by the mm source 
mml or another source nearby. With the current spatial 
resolution it is likely that the outflows South (A) and (C) 
emanate from the same mm continuum condensation. At 
a distance of 2.1 kpc, the beam-size of 3.9" x 3.6" corre- 
sponds to a linear resolution of approximately 7900 AU. 
Thus, it is well possible that a multi-component system 
is evolving within mml and the outflows emanate from 
different membe rs of the system s imilar to observations 
in, e.g., HH288 llCueth et all 1200 ll) or IRAS 05358+3543 
l|Beuther et"aTll20024 ~ 

South (D): The fourth southern outflow labeled (D) 
is further to the west and runs roughly north-south. It 
would already have been possible to identify this flow in 
the single-dish data, but this structure does not stand out 
prominently there (Fig. We do not find a mm contin- 
uum or X-ray source at the center of this outflow (Figs. [21 
& 0) ■ The driving source can well be hidden in the noise 
because the continuum mass sensitivity is > 20 M Q f t]3.1|) . 
The H 2 feature (3) is in the vicinity of this outflow but we 
cannot unambiguously attribute it to South (C) or (D). 

3.2.2. The northern region 

The northern region has a similarly complex morphologi- 
cal structure and we suggest that at least three outflows 
are present there. 

North (E): As for the outflows South (A) and (B), the 
northern outflows (E) and (F) could also be parts of only 
one outflow system, but for those structures, especially 
the CO emission indicates the existence of two overlap- 
ping outflows. Outflow North (E) runs in a collimated way 
roughly in east-west direction and corresponds to the pre- 
viously detected single-dish outflow in this region. Likely, 
the mm source mm5 is the driving source of this outflow. 
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Source 


R.A. 


Dec 


Flux 


Peak 


M 


Nn 2 




[J2000.0] 


[J2000.0] 


[mjy] 


[mjy/beam] 


[M ] 


[10 24 cm" 2 ] 


mml 


19:43:11.24 


23:44:03.4 


29 


19 


132 


4 


mm2 


19:43:11.17 


23:44:07.4 


39 


19 


177 


4 


mm3 


19:43:10.59 


23:44:01.5 


14 


9 


64 


2 


mm4 


19:43:11.17 


23:43:48.0 


13 


8 


59 


2 


mm5 


19:43:10.68 


23:44:58.4 


10 


8 


45 


2 


mm6 


19:43:10.81 


23:45:07.0 


13 


8 


59 


2 



Table 1. The mm sub-sources: positions, peak fluxes, integrated fluxes, derived core masses and column densities. 



Fig. 2. Channel map of the merged PdBI and 30m CO(1-0) observations. The data are binned into 3km/s channels and 
labeled each at the top-left. The arrows, ellipses, letters and markers outline the outflows and mm sources (triangles) 
as discussed in the text. 



Fig. 3. Top: The left and middle panel show the blue (v[5;18] km/s) and red (v[26;47] km/s) CO emission separately 
(merged PdBI and 30 m data) . The right panel then presents the combined outflow emission and marks the outflows 
(arrows, ellipses and letters), mm sources (triangles) and X-ray sources (stars) as discussed in the text. Bottom: The 
left and middle panels show the K' and H2 images from Calar Alto, respectively. The labels in the middle panel mark 
the H2 features discussed in the text. The contours in the right panel again show the PdBI CO emission, this time as 
an overlay on the H2 map. The markers and errors are the same as in the top panel. 



North (F): While the single-dish data suggested a sim- 
ple morphology (Fig. [Q, the red emission between the 
northern and southern core stands out far more promi- 
nently in the PdBI data. Therefore, we suggest that the 
red emission outlines the ellipsoidal walls of a cavity-like 
outflow, as outlined in Figures El & El The blue-shifted 
emission to the north-east shows a cavity-like structure 
as well and supports this interpretation. In Figure the 
cavity structures are best observed at velocities 18, 30, 
33 km/s. In this picture, the H 2 features (4), (5) and (6) 
indicate shocked gas at the interface between the cavities 
of the outflows (B) and (F). There is also a faint H2 emis- 
sion patch (12) in the north which could be part of the 
northern cavity and which is just outside our observed 
PdBI mosaic. The morphology suggests that the driving 
source of outflow (F) is located within the mm source mm5 
but it could also be driven by mm6. 

North (G): The finger- like blue filament pointing north 
has a a red counterpart pointing south, both centering 
around the mm source mm6. This morphology suggests a 
seventh outflow (G) running north-south. Source mm6 is 
approximately at the center of this structure and the only 
mm source where we find high-velocity gas as well as H2 
emission (10) centered on (Figures |21 & EH> ■ This suggest 
that the H 2 feature is directly produced from the underly- 
ing high-velocity jet. At the northern tip, we find weak H2 
emission (11) which can be associated with the blue CO 
emission. Due to the spatial overlap with other outflows, 
it is not clear how far south this North (G) proceeds and 
which part of the emission belongs to North (G) or rather 
to South (A), (B) and (C). There are more H2 emission 
knots (7, 8, 9) associated with the southern part of this 
structure, but those features are difficult to interpret and 



it is not clear whether they are part of North (G) or belong 
to a completely different structure as outlined below. 



3.2.3. Further tentative outflow structures 

In the interface region between the northern and the 
southern core, we find red/blue CO emission and the H 2 
features (7), (8), and (9). 

While it is possible that (7) and (8) are produced by 
outflow North (G), the elongated H2 emission (7) coincides 
with elongated red CO emission. Additionally, Figure ED 
shows blue CO emission at the H 2 knot (7). These obser- 
vations allow the tentative detection of another outflow 
(H) outlined by the dashed line in Figure EH We do not 
detect a mm continuum source there as the outflow driver, 
but, as already mentioned, the continuum sensitivity cor- 
responds only to a mass sensitivity of > 20 M Q , and a 
low-mass source can easily be hidden there. 

In addition to the red CO emission corresponding to 
outflow North (F), there is also strong blue CO emission at 
the interface of the northern and the southern star-forming 
cores. It is not entirely obvious where this structure be- 
longs to, but its north-east to south-west elongation, the 
associated red CO emission and the H 2 knot (9) suggest 
that even a ninth outflow (J) might be depicted by these 
structures. There is no mm source detected in that re- 
gion as well. We stress that the structures (H) and (J) 
are even more ambiguous than the other seven outflows 
discussed before. Higher spatial resolution and sensitivity 
are needed to derive a conclusive picture of the different 
outflow phenomena. 
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3.3. Position-velocity diagrams 

In spite of the seemingly chaotic morphology, the data still 
allow a kinematic interpretation. Therefore, we present 
in Fig. 0] position-velocity diagrams for the five non- 
ellipsoidal outflows (A, C, D, E, G) as well as for the 
tentative structures (H) and (J). North (G) is the only 
outflow with high-velocity gas at its center, for the other 
outflows, the velocities increase with distance from the 
outflow center. This is also the case for the position ve- 
locity diagrams corresponding to (H) and (J), which is 
an additional indication that those features are part of 
outflows. The position- velocity morphologies resemble a 
Hubble-law, which can be interpreted due to the combina- 
tion of momentum conservation of the protostellar winds 
and density decrea ses away from t he driving sources of 
the outflows (e.g., IShu etaD Il99lh . Hydrodynamic jet- 
entrainment models al so result in Hubble-like position- 
veloc ity diagrams (e.g.. lSmith et al.lll997tlDownes fc Ravi 
Il999h . The observed position-velocity diagrams resemble 
those known from low-mass star formation. It has to be 
noted that already the position-velocity diagrams derived 
from the single-dish data show similar features. 



3.4. Outflow parameters 

In addition to morphological and kinematic interpreta- 
tions, we estimate masses, energetics and outflow rates of 
the individual sub- flows. Because there is significant over- 
lap between the various outflows, often we cannot distin- 
guish how much emission corresponds to one or another of 
the overlapping outflows. Thus, some emission features are 
used to derive the parameters for different outflows. In the 
following, we prese nt the outflow paramet er for the single- 
dish observations ijBeuther et all l2002b(l . compare them 
with the outflow parameters derived for the same regions 
from the merged PdBI and 30 m data, and finally discuss 
the outflow characteristics of the individual sub-outflow. 
It has to be noted that the single-dish results presented 
here deviate from the parameters presented in Table 2 of 
iBeuther et all (|2002b). because in the former paper we 
erroneously derived the mass in the red outflow wing of 
IRAS 19410+2336 (south) using a wrong source size. This 
error occure d only in this case, and the other parameters 
of Table 2 in lBeuther et all (|2002bj) are correct. 

Opacity-corrected H 2 column densities and N T in 
both outflow lobes can be calculated by assuming a con- 
stant 13 CO/ 12 CO 1 — line win g ratio throughout t he out- 
flows llOahrit fc Bertoutl fToOfTl . lOhoi et al J jlflfll found 



South (A) 



outh (C) 



an average 13 CO/ 12 CO 2 — 1 line wing ratio around 
0.1 in 7 massive star-forming regions, , corresponding 
to a t( 13 CO 2-1) = 0.1. We adopt this value for 
IRAS 19410+2336 as well, and we assume 30 K as aver- 
age temperature in the outflow. The mass of entrained 
gas M out , the dynamical timescale t and the outflow rate 
-^outflow are calculated via: 



A/out = (Mi x sizc b + N r x size r ) m H2 



—10 - 




10 



15 20 25 30 
Velocity [km/s] 
South (D) 



15 20 25 
Velocity [km/s] 
North (E) 



-10 - 




10 15 20 25 30 
Velocity [km/s] 

(J) 

-40 I-' ' 




15 20 25 30 10 

Velocity [km/s] 



15 20 25 30 
Velocity [km/s] 



Fig. 4. Position-velocity diagrams for the five non- 
ellipsoidal outflows plus the two tentative ones (H) and 
(J). The horizontal line in each diagram corresponds to 
the anticipated outflow center. The latter are only assump- 
tions for (H) and (J). 



M n 



J/2 



with sizcb and size r the sizes of the outflow lobes, mjj 2 



the mass of the H2 molecule, and v n 



and v-n 



the 



maximum velocities observed in each line wing. A more 
detailed description how the outflow parameters are de- 
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termined is given i n IBeuther et alJ l|2002b|) . According to 
ICabrit fc Bertoutl l)l990|) derived masses are accurate to a 
factor 2 to 4, whereas the accuracy of dynamical parame- 
ters are lower, at about a factor 10. 

The first two sections of Table [5] correspond to the 
derived total outflow masses for the single-dish data and 
the merged PdBI data, respectively. The chosen outflow 
region is based on the single-dish data. Comparing the 
masses of entrained gas M ou t from the single-dish and the 
PdBI observations, we find that the merged interferomet- 
ric data compare reasonably well with the single-dish re- 
sults. This approximate agreement is further support that 
the merging of the interferometric and single-dish data, 
and also their relative weighting, has worked properly. 
Additionally, it indicates that single-dish observations are 
adequate to derive global outflow properties. 

More interesting are the parameters of all the sub- 
outflows derived separately. We note that summing up 
the values from all southern and northern outflows does 
not result in the total mass derived before, because the 
morphologies vary too much: e.g., some emission features 
believed from the single-dish data to be part of the south- 
ern flows belong in fact to different flows of the northern 
region. Furthermore, other emission patches are used to 
derive the parameters for different outflows because we 
cannot always distinguish unambiguously where they pri- 
marily belong to. 

The derived outflow masses M out are between I M 
for the smallest outflow South (C) and 12 M Q for North 
(F). The outflow rates are high, of the order a few 
times 10 -4 M Q yr _1 . Assuming momentum driven out- 
flows, a velocity ratio between the entrained gas and the 
underlying jet of approximately 1/20 and a ratio be- 
tween the jet-mass-loss rate and the accretion rate of 
0.3, such outflow rates result in accretion rate estimates 
of the or der 10~ 4 Mm vr -1 (fo r details on the assump- 
tions see IBeuther et all l2002b|) . Accretion rates of that 
order are necessary to overcome the radiation pressure 
and for m massive stars via disk-mediated accretion pro- 
cesses ijStahler et all l200fjjl . The accuracy of mass, age 
and other parameters is not good enough to compare in- 
dividual properties between the different flows in a mean- 
ingful way. 

4. Conclusions 

The new high-spatial resolution CO and H2 observations 
of the high-mass star-forming region IRAS 19410+2336 re- 
solve at least seven (maybe even nine) bipolar outflows. 
The large number of outflows and their interactions are 
an excellent example for the complexity in high-mass star 
formation. The data show that the single-dish observa- 
tions simplify the morphological picture of the region, 
but nevertheless give a good global picture of the region. 
Unfortunately, we also find that even the present interfer- 
ometric spatial resolution is not sufficient to clearly disen- 
tangle the multiple overlapping outflows, different inter- 
pretations for some of the features are well possible. CO 



and H2 emission mostly trace different parts of the same 
outflows, but it is not always possible to attribute all ob- 
servational features unambiguously within the same inter- 
pretation. We note that here we present the - from our 
point of view - most plausible interpretation of the outflow 
morphologies, but it is well possible that follow-up obser- 
vations will result in partly different outflow identifications 
because of the complexity of the region. Nevertheless, the 
whole picture is very suggestive of a scenario, where the 
global outflow properties of massive star-forming regions 
are a superposition of many different outflows emanating 
from the evolving massive cluster. 

The masses and energetics of each sub-flow are, as 
expected, below the values derived from the single-dish 
data. Previous studies of other massive molecular out- 
flows have also shown multiple outflows but the re- 
gions were dominated energeticall y by the most massive 
outflow (e.g., IRAS 05358+3543. IBeuther et all 120023 
IRAS 23033+5951, Wyrowski et al. in prep.). This is less 
clear in the case of IRAS 19410+2336: the northern core is 
dominated by outflow (F), but it does not seem to be the 
case as much in the south, (A) contributes only approxi- 
mately 30% of the mass estimated for the whole southern 
core. The difference is lower if outflows (A) and (B) were 
not two, but are different parts of only one flow. However, 
this question cannot be resolved by the present observa- 
tions. Combining the observations from the literature with 
these new data, we find that global core statistics based on 
single-dish observations are often valid, but the averaging 
aspect of the single-dish data has to be kept in mind. 

Accretion rate estimates based on the outflow energet- 
ics are in the 10 M Q yr _1 regime. Accretion rates of that 
order are important to overcome the radiation pressure of 
the forming massive stars, and thus build up the most 
massive objects via a continuous disk-accretion process. 

Of additional interest are the position-velocity dia- 
grams, which resemble a Hubble-like distri bution as of- 
ten reported in low-mass sta r formation (e.g.. lSmith et al] 
Il997t Downes fc Ravi Il999j) . This result is different to 
findings based on C O single-dish observations from 
Rid ge MooTel l|200lh who report that they do not find 
Hubble-like position- velocity dependence in any of their 11 
studied massive outflow sources. In contrast, resolving the 
two single-dish outflows of IRAS 19410+2336 into at least 
seven outflows with the PdBI indicates that the kinemat- 
ics of the individual outflows resemble the features known 
from low-mass star formation. Furthermore, the H2 emis- 
sion associated with South (C) suggests the presence of 
outflows in a well-collimated jet-like manner, as is found 
in low-mass protostars. Together, this is further support 
to the hypothesis that the physical processes in high-mass 
star formation are qualitatively similar to those found in 
their low-mass counterparts, the main difference being the 
clustered mode of formation and higher accretion and out- 
flow rates (e.g.. iMcKee fc TanlEooa iMaeder k. Behrendl 
20021 IBeuther et al.ll2002blldh . 

Concerning the dust continuum emission, in our pre- 
vious single-dish 1.2 mm observations with a spatial res- 
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Table 2. Outflow results: masses M b (blue), M r (red) and M out (M out = M h + M r ) [M©], momentum p [M© 
km/s], energy E [10 46 erg], size [pc], time t [10 4 yr], mass entrainment rate M out [10 _4 M©/yr], mechanical force F m 
[10 _3 M©km/s/yr] and mechanical luminsoity L m [L©] 



source 


A/ b 


M r Mout p E size t 


Afout 


F m 


L m 






Single-dish observations with the IRAM 30 m in 


CO(2-l) 






South 
North 


11.2 
11.2 


19.8 31.0 683 1.510 47 0.4 19510 
13.4 24.6 541 1.2 10 47 0.5 20420 


1.6 1Q~ 6 
1.2 10" 3 


3.5 10 -2 
2.710 -2 


63 
49 



Merged PdBI and single-dish data of the single-dish regions in CO (1-0) 
South 8?2 WA 18~6 651" 2.6 10 4Y OA 12790 1.5 10~ 3 5.1 10 -2 167 

North 9.3 14.8 24.1 880 3.6 10 47 0.5 13380 1.8 10~ 3 6.6 10~ 2 222 

Merged PdBI and single-dish data of the individual outflows in CO (1-0) 



South (A) 


2.6 


3.0 


5.6 


192 


7.6 10 40 


0.3 


8330 


6.710" 


-4 


2.310" 


-2 


75 


South (B) 


1.7 


1.6 


3.3 


107 


4.1 10 46 


0.4 


11300 


2.910" 


-4 


9.510" 


3 


30 


South (C) 


0.6 


0.7 


1.3 


46 


1.8 10 46 


0.2 


5060 


2.610 


-4 


9.010" 


3 


29 


South (D) 


1.6 


2.6 


4.2 


153 


6.3 10 46 


0.2 


5950 


7.110" 


-4 


2.610" 


-2 


87 


North (E) 


3.5 


4.2 


7.7 


265 


1.1 10 47 


0.4 


11300 


6.810" 


-4 


2.310" 


-2 


76 


North (F) 


4.5 


8.1 


12.6 


472 


2.0 10 47 


0.6 


17840 


7.110" 


-4 


2.610" 


-2 


91 


North (G) 


1.9 


6.8 


8.7 


359 


1.6 10 47 


0.6 


16950 


5.1 10" 


-4 


2.110" 


-2 


77 



olution of 11", we found that the power-law density dis- 
tribution does not continue to the very center, and we 
derived an inner breakpoint at 5" for the southern core 
from where on we mod eled the density distribution as flat 
l)Beuther et all l2002afk We suggested that this flattening 
is due to fragmentation of the large-scale core into many 
sub-sources. The new interfcromctric data directly confirm 
this hypothesis, and we observe the fragmentation into 
sub-cores (Fig. [I). The spatial separation of, e.g., mml 
and mm2 is ~ 5". 

To summarize, the presented observations show the 
complex interaction of several molecular outflows within 
the massive star-forming region IRAS 19410+2336. In 
spite of the complexity, the analysis of the data supports 
the idea of massive star formation taking place in a clus- 
tered mode by similar physical processes as known for 
low-mass star formation, involving accretion through disks 
producing collimated outflows. 

Finally, we note that follow-up observations of this re- 
gion are necessary to understand the interplay of many 
molecular outflows within one star-forming region more 
properly. In spite of the similarity to low-mass star forma- 
tion, the feed-back processes between many protostars in a 
high-density environment should also produce significant 
differences, IRAS 19410+2336 can be used as a template 
region for such studies. Regarding molecular outflows, dif- 
ferent molecular observations are needed, especially inter- 
esting is SiO which is known to be a shock and jet tracer. 
Furthermore, we are still resolution limited, and it will be 
crucial to re-observe the region in CO(2-l) with signifi- 
cantly higher spatial resolution. 
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